The conductance of superconducting-normal hybrid structures 
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The dc conductance of normal-superconducting hybrid structures is discussed. It is shown that 
since the Bogoliubov-DeGennes (BDG) equation does not conserve charge, its application to create 
a Landauer-type approach for the conductance of the NSN system is problematic. We 'mend' this 
deficiency by calculating the conductance from the Kubo formula for a ring configuration (for this 
geometry the solutions of the BDG equation conserve charge). We show that the presence of a 
superconductor segment within an otherwise normal metal may reduce the overall conductance of 
the composite structure. This reduction enhances the tendency of the NS composite to become 
insulating. 
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I. INTRODUCTION 

In a seminal paper published a long time ago, Blonder, 
Tinkham, and Klapwijki (BTK) calculated the conduc- 
tance of a normal (N)-superconducting (S) interface as 
a function of the interface transparency. In particular 
they showed that at zero temperature (T = 0) and for 
electrons at the Fermi energy that conductance is given 
by 
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where jS'^'^p is the "normal" reflection coefficient of the 
interface, while |S''''^p is the reflection coefficient for the 
Andreev processes^ which reflect electron-like excita- 
tions as hole-like ones. For simplicity, both the nor- 
mal and the superconductor regions were taken to be 
free of any impurity scattering (except at the interface) . 
Then, when the interface is perfectly transparent, i.e., 
1 5'='= 1 2 = and 1 5''^ 1 2 = 1, the value of G^g is t^ice 
that of the "normal" quantum limit of the conductance. 
However, when a large enough barrier exists at the in- 
terface between the superconducting and the normal re- 
gions, Gj^g becomes much smaller than the conductance 
obtained when the superconductor is made normal. This 
is the simplest example where superconductivity in a 
part of a system reduces its overall conductance. (The 
more complicated many-channel, disordered case is not 
addressed in this paper.) In this paper we consider the 
more subtle NSN combination, and demonstrate a simi- 
lar effect: with a large enough barrier at even one of the 
NS interfaces, the appearance of superconductivity in the 
S region reduces the overall conductance. 

Some of the motivation for the present work comes 
from our wish to understand why the mixed NS bulk 
composite structure is often insulating at T = and the 
superconducting phase of a thin film goes over (with in- 
creasing disorder or decreasing film thickness) directly 



into the insulating rather than to the normal-conducting 
phase-^;-i It is interesting that often the activation energy 
for the conductance of the insulating phase is given by the 
superconductor gap of the superconducting component 4 
Thus, the superconducting component plays the role of 
an additional barrier between the normal segments! The 
charge-vortex duality^ for a system of charged bosons of 
course explains the insulating phase as dual to the super- 
conducting one. Vortex localization yields zero resistance 
and charge localization yields zero conductance at T = 0. 
Our purpose is to provide a heuristic understanding of 
how the charge localization is established. Reducing the 
small-scale conductance of the system pushes it towards 
the insulating state. The NSN system is the simplest mi- 
croscopic element of the NS network. We find that it al- 
ready presents nontrivial theoretical questions having to 
do with a seeming deficiency of the Bogoliubov-DeGennes 
(BDG) formulation. 

Blonder et alX employed the BDG equation for the 
quasiparticle excitations in the superconducting region, 
assuming that the energy gap A which vanishes in the 
normal part does not vary spatially in the superconduc- 
tor. Since in most situations the superconducting coher- 
ence length ^ is much larger than the Fermi wavelength 
and much smaller than the length of the S-rcgion, this 
assumption seems quite harmless. Consequently, it is 
widely accepted that the use of the BDG equation, with- 
out attempting to compute the (complex, in general) su- 
perconducting order-parameter self-consistently, is valid 
for many hybrid structures. 

A particularly important issue emphasized by BTK 
concerns the conversion of the normal current into 
a super-current at the NS interface, ensuring charge- 
current conservation over the entire structure (see also 
Ref. 6). Blonder et al} showed that the normal charge- 
current (which is distinct from the quasiparticle current) 
entering the superconductor decays, but concomitantly 
a super-current grows up gradually, until very far inside 



S the normal current disappears completely, and the en- 
tire charge is carried away by Cooper pairs. One might 
wonder what happens to this scenario when the super- 
conductor has a finite width and is not infinite as in the 
BTK case. It turns out that within the BDG formulation 
the super-current in the S region does not properly con- 
vert back to normal current in the second N region. We 
are not aware of a way to correct for this deficiency (leav- 
ing aside the possibility mentioned above to compute the 
order-parameter self-consistently) . In this paper we will 
circumvent this problem by using a particular geometry. 
Not surprisingly, transport through hybrid NS struc- 
tures has been addressed before, even prior^ to the pub- 
lication of the paper by BTK. Indeed, as Eq. ^ bears a 
strong resemblance to the Landauer formula* for coher- 
ent transport, several modern treatments (repre sentative 
examples are to be found in Refs. IS flOlllllll^ ) employ 
scattering theory within the Landauer picture in an at- 
tempt to extend the BTK result ([T|) to more than a sin- 
gle NS interface. The scattering formalism is also used 
to study the current fluctuations in mesoscopic systems 
with Andreev reflections J^i-^^ However, there is a vicious 
caveat in this approach: the BDG equation, while con- 
serving the number of quasiparticles, does not conserve 
charge. The BDG formulation follows, in this respect, 
the similar deficiency of the BCS approach upon which 
the BDG equation is based. This does not cause any 
harm when the formulation is used for a bulk supercon- 
ductor, as had been done by BTK (see also Ref. il3), or 
when the system has the shape of a closed ring, where 
the boundary conditions enforce current conservation.^^ 



However, when the size of the superconducting segment 
is finite, the non-conservation of the charge leads to cur- 
rent vs. chemical-potential-difference relations (within 
the linear response regime) that depend on the chemical 
potential of the superconductor, as opposed to the situ- 
ation in the normal phasci^ii^ Indeed, the conductances 
of various hybrid structures have been determined within 
scattering theory (when the superconducting is "float- 
ing" ) by fixing the chemical potential on S such that cur- 
rent is conservedi^ii^ One would have assumed naively 
that this procedure cures the problem mentioned above. 
It turns out, however, that it does not: as discussed by 
Anantram and Dattapii (see also Ref. [l3 ). applying the 
very same procedure to the calculation of the current fluc- 
tuations (the power spectrum of the noise) violates the 
Johnson-Nyquist relationships. We show in Sec. |V] that 
this 'floating-potential' approach also does not produce 
the doubling of the conductance [see Eq. ^ above] when 
the NS interface becomes transparent. It approaches, 
however, the conductance found by the Kubo formula in 
the limit where the barrier is strong and its transparency 
is low. 

It follows that in order to use the BDG equation for 
the calculation of the conductance of a mesoscopic system 
containing a finite-length superconducting segment, one 
needs charge-conserving solutions of that equation. Such 
solutions arise naturally when the NSN structure (the 
left panel in Fig. [T]) is closed to form a "ring" (the right 
panel in Fig. [1]). One can use those solutions in the Kubo 
formula for the conductivity of a large system. This is 
the approach adopted in this paper. 
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FIG. 1: Left panel: tlie 'open' NSN structure, wliere the perpendicular lines represent potential barriers (the explicit model 
treated in Sec. IIVI assumes that the right interface is clean). Right panel: the corresponding ring NSN system. The incoming 
and outgoing amplitudes of the electron-like and the hole-like waves on both sides of the superconductor are marked by arrows. 
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We address the simplest problem of a one-dimensional 
NSN structure at zero temperature (we do not discuss 



here the complications arising in the multi-channel case 
including disorder) . After discussing the conductance be- 



tween the two N segments (left panel in Fig. [T]) , and deal- 
ing with the nontrivial problem arising from the above- 
mentioned deficiency of the BDG picture, we find that 
the NSN conductance can indeed decrease when, for ex- 
ample, the superconducting component becomes longer. 
We are not able to treat here the behavior of the two- 
dimensional or the three-dimensional NSN arrays. Suf- 
fice it to say that when the "small scale" resistance of the 
elementary building block increases, the tendency for lo- 
calization at larger scales becomes stronger. 



The middle conductor, of a finite length, is kept floating, 
(i.e., it is not connected to any reservoir), and will ac- 
quire a chemical potential /i„. Clearly, the right-going 
currents to the left of the middle segment, II, and to its 
right, Iji, are given by^^ 

2e / 

^L = -j- ( (1 - 'R-){^J■L - M«) - T'ifiR - ^n) 

/fl = ^ (-(1 - n'){^iR - M„) + T{^iL - Mn) ) ■ (3) 



II. DIFFICULTIES WITH THE 

LANDAUER-TYPE FORMULATION FOR THE 

NSN CONDUCTANCE 

The subtleties involved in producing a consistent 
Landauer-type formula for the conductance of an NSN 
structure within the BDG formulation are best explained 
by considering the simplest single-mode, two-terminal 
conflguration at zero temperature. In other words, we 
assume that there is no scattering in the system except 
for the potential barriers at the interfaces, such that the 
transverse channel modes are not mixed (and their in- 
dices can be omitted). 

We start with the purely normal case, as shown in the 
left panel of Fig. [1] except that the S section is replaced 
by a normal one (for example, by letting its gap approach 
zero) . In this case there is no need to treat electrons and 
holes concomitantly, as there are no Andreev processes. 
It is enough to use only electrons (or only holes). We de- 
note the reflection coefficient for an electron coming from 
the left by 7^, and the one for an electron coming from 
the right by TV . Likewise, the transmission coefficient 
from the left to the right is denoted by T, and the one 
from the right to the left is T'. Unitarity (particle con- 
servation, which is also charge conservation in this case) 
implies 



7^ + r = 7^' + r' = 1 



(2) 



Time-reversal symmetry implies further that T — T' , 
and hence TZ = TZ' . Next we assign to the left conductor 
a chemical potential fiL and to the right one a chemical 
potential ^r. In the linear response regime, /i^, — /i_r -^ 0. 



From the unitarity condition ([2]) and time-reversal sym- 
metry, it follows that II = Ir = I = i'2e/h)T(nL — Mij), 
independently of the value of //„ (which, in fact, drops out 
of the two equations). The well-known Landauer formula 
for the conductance, 



G = {2e^/h)T 



(4) 



is immediately obtained, and is independent of /z„, as it 
should be.— 

When the middle section is a superconductor, further 
Andreev-type processes become possible. An electron 
can be refiected/transmitted as a hole, and vice versa. 
For an electron incident from the left, the probabilities 
for the Andreev refiection and transmission processes are 
denoted TZx and 7a, respectively. The corresponding 
quantities for an electron coming from the right are TZ'p^ 



and 7](. The unitary condition ([2]) is now replaced by 



7^ + r + 7^A + Ta = 7^' + r' + 7^'^ 



rk 



(5) 



However, the charge conservation condition now reads 



,17 



7^ + r - 7^A - Ta = 7^' + r' - 7^'A - t'a = 1 . (6) 

The two conditions, Eqs. ([5]) and ([6]), are not compatible 
whenever the Andreev probabilities do not vanish, except 
for a ring geometry (where their consistency is enforced). 
Moreover, while Eq. (O always holds for the solutions of 
the BDG equation, Eq. ([6]) does not! 

The expressions for the currents [see Eqs. ([3])] now 
become (note that group velocity of the holes is opposite 
to that of the electrons) 
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((1 - 7^ + 7^A)(/iL - /is) - {T' - TDifiR - fis)) , 



iR^=^[-{i-n' + 7^^)(/ifl - us) + (T - TA){fiL - Ms) ) , 



(7) 
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where ^s is the chemical potential on the superconduc- 
tor. We note that because charge conservation does not 
hold [see Eq. ©], Hs does not drop out of these equa- 



tions. Its value is relevant. Equations ([7]) are of the same 
form as Eqs. (2) of Takanc and Ebisawai^ these authors 
determine ns so that II = Ir- It is then possible to 



obtain a conductance from the current-to voltage ratio, 
as was done in Ref. l9|. We reproduce their resuh for the 
NSN conductance in Sec. |Vl see Eq. ([57| there. It is 
disturbing, however, that the determined value of iis is 
relevant. Moreover, for the result for the conductance 
[Eq. ([57|) ] to be fully satisfactory, it should agree with 
the linear response, Kubo, formula for the related large 
ring geometry. This is the case for the normal Landauer 
formula, but, in general not for Eq. (I57p (see Sec. IVlfor 
more details). 

The above considerations, in particular, Eqs. ([5]) and 
([S]), can be put on a more general basis. By imposing 
the appropriate boundary conditions on the plane-wave 
solutions of the BDG equation it is possible to derive the 
scattering matrix, S, of the NSN structure. This (4x4) 
matrix relates the amplitudes of the incoming waves to 
those of the outgoing ones, (see left panel in Fig. [T]) 
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Herefi^ the incoming amplitudes are 

C;,,= [ct{N),c-{N'),c^{N),c+{N') 
and the outgoing ones are 



(8) 



(9) 



In Eqs. dH) and (fTO| . c^ ,^{N) denotes the amplitude 
of an electron-like (hole-like) excitation with a positive 
wave vector k^^h incident from the left normal side while 
c~ f^{N) refers to the waves having negative wave vectors. 
Since the BDG equation conserves the number of quasi- 
particles, the scattering matrix S is necessarily unitary, 
and therefore 



''out ''out ''in ''in 



(11) 



However, conservation of the charge currentiiii requires 
[seeEq. ©] 

\ctiN)f-\c-iN)f + \ctiN)f~\ci^{N)\' 

= \ct{N')\' - \c-{N')\' + \4{N')\^ - |c^(7V')P • 

(12) 

Comparing Eqs. pT]) and IT^ . we see that they imply 

|c+(7V)P + |c-(iV')P = \ct{N')\' + \c-{N)\^ , 
\cJ.iN)\' + \ctiN')\' - |c+(iV)p + |c,;(iV')P . (13) 

Namely, the sum of the amplitudes squared of the incom- 
ing electron-like excitations is equal to the sum of the am- 
plitudes squared of the outgoing electron-like excitations, 
and so is the situation for the hole-like ones. These condi- 
tions are the same as those that would have been derived 
from Eqs. ^ and ^, had we required that both condi- 
tions should be satisfied together. This always holds for 
a normal system, in which these two types of quasipar- 
ticles are not mixed. However, in a superconductor the 



Andreev processes mix the hole-like with the electron- 
like excitations, thus violating the conditions p^ for 
general hybrid structures with a finite-size S segment. 
An exception is the ring geometry. There, (see Fig. [T|) 
the ratios c^ h^N) / c^ ^^[N') are necessarily phase factors, 
and so are the ratios cJ^(iV)/c^^(iV'). As a result, the 
plane-wave solutions of the BDG equations for the ring 
geometry do satisfy both conditions P^ . namely, these 
solutions conserve charge. One may therefore employ 
the BDG equation for the ring geometry in the context 
of the Kubo formulation to calculate the conductance of 
the NSN structure. 



III. THE KUBO FORMULA FOR A LARGE 
RING 

For an infinite system, the Kubo-type conductivity at 
frequency uj may be most easily obtained by calculat- 
ing, using the golden rule, the power absorbed by the 
system from a classical monochromatic electromagnetic 
field. We consider for simplicity noninteracting fermions 
(or Fermi quasiparticles), and focus on the Uxx compo- 
nent of the conductivity, 
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(14) 



Here, |j) and \£) are the quasi-electron states and /(e^) 
and f{ti) their populations. In Eq. p^ . V is the volume 
of the system, to be sent to infinity at the end of the cal- 
culation, at which stage the summations over the states 
are replaced by integrations with the densities of states. 
The a;— component of the velocity operator is denoted Vx- 
For the case of a normal (i.e., non superconducting) scat- 
terer with infinite leads the equivalence of the Kubo and 
the two-terminal Landauer approaches has been estab- 
lished in Refs. [H andHl. 

The assumption of an infinite system is crucial in or- 
der to have a continuum of states. An isolated finite 
system with a truly discrete spectrum does not in fact 
absorb energy from the monochromatic field. In order 
to obtain a finite conductivity for a finite large system, 
it has to be (and it is, in most real situations) coupled 
to a very large heat bath. For example, to an assembly 
of thermal phonons. This enables energy to be trans- 
ferred from the electromagnetic field into the bath via 
the small electronic system. For a weak enough interac- 
tion with the bath, one may say that the discrete levels of 
the system have acquired finite widths, ?/j. It then makes 
sense to write down Eq. (J14p with the levels having a fi- 
nite width (or with an imaginary part to the frequency 
w, which will amount to a non- monochromatic driving 
field). This procedure has been discussed, including the 
dc limit (Rew ^- 0), by Thouless and Kirkpatrick,^° fol- 
lowin g C zychoU and Kramer^ and used for example in 
Ref. [23, see also Ref. il9. It is postulated, and can 



be demonstrated in typical cases, that once the rjj's are 
larger than the level spacing near the Fermi energy, but 
much smaller than all other relevant energy scales in the 
problem, this procedure yields the physically relevant 
low-frequency conductance of the system. 

It hence follows that the a; ^ conductance is ob- 
tained upon transforming the summations in Eq. (J14p 
into energy integrations. This allows one to approxi- 
mate /(cf -I- hu) — f{ee) ~ hijf'{ee). Focusing on our 
one-dimensional configuration, we take x along the ring 
circumference, and replace the volume of the system by 
its length, d. Since the (one-dimensional) conductance is 
related to the conductivity by G = cr/d, we recover the 
Kubo-Greenwood-type formula at low frequencies 
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(15) 
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where the sum is over the (almost) degenerate initial 
states and over the (almost) degenerate final states, 
within the narrow range hio (— > 0) above those initial 
states, and all states are at about the Fermi energy. One 
might also say that the cancellation of the frequency [see 
Eq. P^ ] is caused by the fact that the initial state was, 
at T = 0, within hu of the Fermi level. In Eq. P^ . 
the matrix element squared of the velocity was replaced 
by its typical value in the small relevant energy window 
around the Fermi energy. The double sum of Eq. (J14p 
gave rise to two factors of the single-particle density of 
states (per unit energy, per unit length, and per spin), i/. 



u = 1/(/iuf) , 



(16) 



in the one-dimensional system. Comparing Eq. (jisp with 
the "traditional" Landauer formula ([4]), we find that in 
the Kubo approach the total transmission is replaced by 
the appropriate sum over the velocity matrix elements 
squared, i.e.. 
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(17) 



It is instructive to review the way the Kubo formula in 
its form (fT7|) produces the Landauer result for the usual 
two-probe geometry. We consider initial left-going scat- 
tering states. These are degenerate with the right-going 
ones. This degeneracy gives a factor of two in the final 
result, to which the spin degeneracy adds another fac- 
tor of two. Each such state will have a matrix element, 
(1 — 7?. + 7")/2 = T, with the appropriate final left-going 
scattering state, and rt with the final right-going scat- 
tering state (r and t are the reflection and transmission 
amplitudes, respectively). Adding the absolute values 
squared together yields T (note the cancellation of the T^ 
term!). Introducing the above degeneracy factors gives 



deg I 



1 1 = ivpT, and thus reproduces the Landauer 
result, Eq. ([4]) above. 

For the ring geometry the states are stationary and 
normalizable. Taking as a representative example a nor- 
mal ring with a single delta-function potential, one finds 



that the ratio of the amplitudes of the clockwise mov- 
ing wave and anticlockwise moving one is a phase factor, 
exp[i(/)e], where on the Fermi energy exp[i(/)e] = —1 or 
(1 — iC)/(l + *C)- Here ( is the strength of the delta 
function potential, with the corresponding transmission 
T = 1/(1 + C^). The velocity matrix elements are then 
vf/(1 ± iC)i E^nd thus together with the spin degener- 
acy reproduce the Landauer formula ^. We give more 
details in the next section, which is devoted to the eval- 
uation of the states and the current matrix elements for 
an NS ring, and the case of an entirely normal ring is 
treated as a limiting case. 



IV. THE VELOCITY MATRIX ELEMENTS 

Here we compute the matrix elements of the velocity 
operator, which are used in the Kubo formula (J17p for the 
conductance. We follow BTK in assuming that the entire 
scattering takes place only at the NS interfaces, and that 
the pair-potential A is finite and spatially-invariant in 
the superconducting region, and vanishes in the normal 
one. Since then there is no channel mixing, the problem 
becomes effectively one-dimensional, and the quasipar- 
ticles are described by the one-dimensional Bogoliubov- 
DeGennes equation, (we use in this section units in which 
h=l) 
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*(j:) = e^(a;) . (18) 



Note that this equation takes into account the two possi- 
ble spin directions pertaining to a certain energy e (mea- 
sured from the Fermi level). In N, where A = 0, the 
solutions of Eq. p8|) are 
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(19) 



with the wave vectors 



k^^f, = ^/2m{EF ± e) ~ fcp ± e/wp • (20) 

In the superconducting segment the solutions are 
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Here, 
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and 
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Up 



Vt= V^2I:a2 , e> A , 
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(21) 
(22) 

(23) 
(24) 



The factor [e/l^]^/^ compensates for the different group 
velocity of the quasiparticles in the superconductor 
[de/dq — {q/m)[n/e)\, and it multiphes the usual coher- 
ence factors u and v, u — [e/fij^/^u and v — [e/n,]^^'^v. 
The amplitudes cff^{N) (see Fig. [T]) are the coefficients 
of the waves exp[±iA:e,?i2;]. Analogous amplitudes are de- 
fined for the waves in the superconducting segment. 
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For simplicity, we assume that the (left) NS interface at 
X = (see Fig. [T]) is represented by a delta-function po- 
tential, XS{x), of strength C = X/vp. Then the boundary 
conditions^'' arc the continuity of the wave functions and 
the discontinuity (of magnitude C.) of their derivatives, 
leading to the following relations among the amplitudes 
of the N region and those of the S one. 
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The other NS interface, located at x = ds, is assumed to be perfectly transparent, and then the boundary conditions 
are the continuity of the wave functions and their derivatives. When the system has the shape of a ring, in which the 
length of the normal segment is cIn, these boundary conditions are 
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where 
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is 



JQds/vp 



(27) 



Without loss of generality we may choose exp{ikpd) = 1, 
where d = d^ + ds is the total length of the ring. Then 
fcp disappears from the boundary conditions. 

Upon eliminating the S-region amplitudes, one obtains 
the equation which determines the allowed eigenenergies 
of the ring, and the ratios among the amplitudes of the 
normal region for each such energy. 
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Here, 



X = 7jyMcos(r2ds/wF) - i-^ sm^nds/vp) 
Y — 2ijpjuv sm{flds /vp) , 



(29) 



such that |Xp- |y|2 = 1 for both e > A and e < A. The 
allowed eigenenergies are given by the vanishing of the 
determinant of the matrix in Eq. (|^5|) . The zeroes of the 
determinant define the families of possible eigenenergies. 



which are rather dense when the size of the entire system 
is large. When (^ ^ 0, the eigenvectors of the matrix 
([28)1 are such that the ratios of the clockwise electron 
(hole) waves to the anticlockwise electron (hole) ones (see 
Fig. [T]) for each of the families of eigenenergies are phase 
factors, 
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In other words, at any finite value C, of the barrier, there 
is a perfect reflection of the electron and the hole waves 
(in the ring geometry). On the other hand, the ratios of 
the hole amplitudes to the electron ones obey 
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such that the phase of P is {(j)^ — (ph)/^, 
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It is illuminating to consider Eq. (p8|l and its solu- 
tions in the limit of very high energies, e ^ A, where 
the superconducting order parameter A becomes irrele- 
vant, and the entire system behaves as if it were normal. 
Then [see Eqs. dM])] X = exp[ied/vF] and F = 0, and 
Eq. (pS)) separates into two independent blocks, for the 
electron-like excitations, and for the hole-like ones. The 



eigenenergies are given by 



ed/v-p 1 



i + < 



for the electron waves 



zed/i,F ^ ^ ^^ I — < ^^j. ^j^^ ^-^ (.33^ 

1 + < ' ' V ; 

with the corresponding phase ratios 



o«0c 



-1 , or 



l + < 



(34) 



(In this hmit P, the ratio of the hole amplitude to the 
electron amplitude, is not defined.) We show below that 
these are the phase factors exp[i(/)e_;i] which determine the 
conductance @ when calculated from the Kubo formula 

m- 

In the other extreme limit of sub-gap energies, e <C A, 
one approximates^ [see Eq. (l23t ] 



(35) 



rj _ A _ 1 

and consequently [see Eqs. (|29p ] 

X ~ 7j^icosh(d5/0 , ^^ - -*7Arsinh(ds/C) , (36) 

where ^ is the coherence length in the superconductor. 
One then finds a quadratic equation for cos{edtq /v-p). We 
do not present explicit expressions for the solutions and 
the amplitude ratios since they are rather cumbersome. 
The next step in this calculation is to find the normal- 
ization of the wave functions, using 

/ "" dx\^s{.x)\'' + / dx\^N{x)\'' = 1 . (37) 

JO J ds 

In the N region 

[*^(x)[2 = [c+(iV)[2 + [c-(iV)[2 

+ (e2''=^"c-(Ar))*c+(iV) +cc) -t- (e ^ h) . (38) 

When djv is large, such that the oscillatory terms (in 
k^disi) can be ignored, the contribution of the normal 
part to the normalization integral becomes 

rdx\^Nix)\'^dN(\ct{N)\^ + \c-{N)\' 

•Ids ^ 

+ \ci{N)\^ + \cn{N)\') . (39) 



The calculation of the contribution to the normalization 
coming from the S region is more subtle, since the wave 
vectors can have an imaginary part [see Eqs. (I22p and 
Disregarding terms oscillating with kpds, we find 



■ (i"i*-n)x 



(n-ti*)x 



2\„«- 



i\c:{S)\' + \ct{S)ne 



-h {u*v + uv*) I e 



■i{c~{s)rc-{s) 

+ «(^))*4(^)) + cc) . (40) 



This rather complicated result reflects the fact (specif- 
ically, its second part) that in the superconductor the 
electron waves are mixed with the hole ones. However, 
at very large energies, e 3> A, or at very small ones, 
e <C A, the mixing term, {u*v + uv*), vanishes. In the 
following, we confine ourselves to these two limits. In 
the high-energies limit the normalization of either the 
clockwise waves or the anticlockwise ones is simply v2d, 
where d is the total length. In the low-energies limit we 
find 



ds 



dx\^s{x)\' = |sinh(ds/0(|c+(iV)|^ + |c-(iV)[2) 



{e'^^'lMj' 



-ds/i 



M, 



where 



Ma = In - ^InP , 
Mb = 1n^ + ^nP , 



(41) 



(42) 



and P is given by Eq. (PT|) . 

Having fully determined the wave functions, it remains 
to compute the matrix elements of the velocity, 

with the indices j and £ enumerating the various eigen- 
functions. As in the calculation of the normalization, 
here again there are contributions from the normal and 
from the superconducting regions. In each region we dis- 
card the oscillatory terms, those which involve k-pd^ or 
kpds. 

The contribution of the normal part to the integral in 
Eq. (gSl) reads 



..N 



dnk 



ATKp 



((c+(iV))*(c+(iV)), - {c-{N))*{c-{N)), 

+ {cl{m*M{m,-{cl{N))){cl{N))^ . (44) 



In the high-energies limit, e 3> A, the contribution of the 
superconducting segment to the integration is the same 
as p4|) (with the arguments iV replaced by S, and d^ 
replaced by ds)- The contribution of the S part in the 
limit of very low energies, e <C A, is 

X ((c+(7V))*(c+(iV)),M^.,-(c,-(iV))*(cj(iV)),M;^) , 

(45) 

where we have denoted [see Eqs. (|42l) ] 

M^, = e'^-^/«(Mj*(Mj, + e~''-/«(M,)*(A4), . (46) 

It is again useful to examine the limit of high energies, 
e ^ A, where the entire ring behaves as if it were normal. 
Then, the electron- and the hole- like waves are separated. 
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The spectrum and the ampHtude ratios are given by Eqs. 
([55)1 and ([M]) . and the normaUzation for each species is 
vM. The matrix elements of the velocity are simply 



"jt 



¥(4):(4).(i 



), 



(47) 



and an analogous result is obtained for the contribution 
of the hole waves. Obviously, the diagonal ones vanish. 
The non-diagonal ones give [k-p /'m)/{\ ± iQ), and conse- 
quently 



El 

dcg 



|2 = 4v2r , T - 



1 



i + C' 



(48) 



Note that the non-vanishing matrix elements arise from 
the phase factor between waves belonging to the same 
species but moving along opposite directions. Hence, the 
Kubo formulation for the ring geometry reproduces the 
Landauer result for the dc conductance. 

Another illuminating limit is when C, vanishes, and 
both NS interfaces (see Fig. [T]) are perfectly transparent, 
in which case the clockwise and the anticlockwise ampli- 
tudes are independent. The matrix elements of the ve- 
locity, for sub-gap energies, are (for either the clockwise- 
moving or the anticlockwise-moving excitations) 



VF , 



V 



v-p- 



dn 



^' ^' dN + etanh(ds/0 
sinh(ds/0[sinh(ds/C)+«] 



cosh" {ds/O 



(49) 



(It is interesting to note that the off-diagonal matrix el- 
ements are coming from the N region alone.) Hence the 
contribution of both the clockwise waves and the anti- 
clockwise waves is 



deg 



4w^l-f 



dNtai\\i{ds / i) 



dN + ^tanh(d5/^) 



(50) 



Thus, when ds/£, tends to zero (namely, in the absence of 
the superconductor) the result approaches the Landauer 
formula for a transparent barrier, cf. Eq. (|48p . On the 
other hand, when d^ > ds ^ (,, our result (|50p tends to 
the one found by BTK^ (for a clean interface), namely, 
it is twice the value of the quantum conductance. 



Unfortunately, the explicit expressions for the velocity 
matrix elements at low energies for general values of C and 
ds/£, are rather complicated. Consequently, we present 
the results of the calculations only graphically, see Figs. 
[2] and [3] The figures show the conductance [divided by 
2e^//i, see Eq. (|17|) ] as a function of the ratio ds/^ for 
various values of the interface transmission, T = 1/(1 + 
C^), and as a function of that transmission, for various 
values of the size of the superconductor segment, ds/C- 



dpj=d^ 




dN=d.; 



dsl? 




FIG. 2: Left panel: Conductance vs ds/S, for several values of the transmission. Right panel: Conductance vs the transmission 
for several values of ds/C- Here the length of the S region equals that of the N region (djv = ds). All information pertaining 
to values of ds smaller than ^ is presented by dashed curves. 



r 



Figure [5] presents the results for the case dN = ds, i.e., 
the segments N and S of the ring are of equal lengths. 
In the left panel the conductance is plotted as a function 



of ds/C) for values of T ranging between 1 (the upmost 
curve) and 0.06, (the lowest-lying one). The main feature 
of these curves is the variation of their slope as the NS 



barrier becomes less and less transparent. For T = 1 the 
conductance increases with the length of the supercon- 
ductor (until it is double that of the normal system in 
the BTK limit where ds/S, -^ c»). As the transparency 
decreases, the conductance, albeit increasing with ds/S, 
becomes smaller, until at about T ~ 0.8 it changes its 
slope and begins deceasing as the size of the supercon- 
ductor is increased. The same characteristic behavior is 
obtained when the size of the normal part largely exceeds 
that of the superconductor, as is depicted in the left panel 
of Fig. [31 The right panels in both Figs. [2] and [3] show 



the (normalized) conductance as a function of the barrier 
transparency for various values of the superconducting 
size, ds/^, ranging between 0.01 (almost a straight line) 
and 20 (parabolic curve). Here one observes that the con- 
ductance is linear in the barrier transmission as long as 
the superconducting is small enough, and then becomes 
quadratic in T, for large values of ds/S,- It should be 
noted, however, that the use of the BDG approach for 
ds/S, ^ 1 is dubious. For this region, we have presented 
all information pertaining to such values by the dotted 
curves. 



J 



drj=9d^ 



=9d. 




2.0 r 



rfs/f 




FIG. 3: Left panel: Conductance vs ds/£, for several values of the transmission. Right panel: Conductance vs the transmission 
for several values of ds/(,- Here the length of the S region is smaller than the length of the N region (djv/ds ~ 9). All 
information pertaining to values of ds smaller than ^ is presented by dashed curves. 



r 



V. DISCUSSION 

As is described in Sees. UlandHIl several previous cal- 
culations aiming to determine the conductance of hybrid 
normal-superconducting structures are based on the scat- 
tering matrix for the quasiparticles, as derived from the 
BDG equationiSiiiSiii Our reservations regarding this pro- 
cedure are explained in Sec. [TTl Nonetheless, it is inter- 
esting to compare the conductance found from the Kubo 
formula and the one derived after fixing the chemical po- 
tential of the superconductor, as explained in Sec. [TTl 
Here we carry out this comparison for the model system 
of Sec. [IVI 

The scattering matrix of the NSN junction (see the 
left panel in Fig. [1]) is a function of the energy e. For 
our purposes here it suffices to derive it for zero energy. 



i.e. on the Fermi level. This derivation is accomplished 
by eliminating the amplitudes of the waves within the 
superconductor, using the boundary conditions ([^5)) . and 
the boundary conditions at the (clean) interface between 
the superconductor and the second normal layer, denoted 
N' [note that when e = 0, 7Ar = 1, see Eq. ([77)1 ] 



ctiN') 1 




U'Jg "fg^V 




\4{S) 


C7{N') 




"75^ ^gV 




cliS) 


ctXN') 




vis Is^'^ 




4is) 


^1{N')\ 




. «7s^ 7s". 




l^hiS) 



(51) 



As a result, the scattering matrix as defined in Eq. 
takes the form 



J 



r c-{N) 1 




4{N') 


1 


ctiN) 


^ D 


[cj;{N')\ 





*C(i 


-zC)(c2-t- 


s') 


(l-zC)c 




isc 




-Cs 


c 


(1 - <) 




-zC(i - *C) 




(s 




isc{l + 2^ 




isc 




c^ 


<(1 


+ iO{c^ + 


■s') 


c{l + iC) 




-Cs 




iscil + 2C2) 




c(l + lO 




<(! + <) 



4{N) 
c-{N') 



(52) 
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where 

i^ = (l + C')c2 + C's' , (53) 

and in order to shorten the notations we have denoted 



sinh((is/0 , c ee cosh(ds/^) 



(54) 



Referring to the notations introduced in Sec. |TT1 we 
find from Eq. (|52p that for an electron-hke incident from 
the left 



^ ^ c!(i + c!i(£i+£!)! ^ ^^c2(i + n 



D2 



D2 



7^A = 



D2 ' 



Ta 



£)2 ' 



(55) 



while for an electron-like wave coining from the right the 
corresponding probabilities are 



7^' 



n\ = 



(2(1 + c^) 



D2 

S2c2(l- 



r' 



\^+e) 



2C 



2\2 



D2 



n = 



D2 



£12 



(56) 



It is easy to verify that the conditions for quasiparticle- 
number conservation, Eq. ([5]), are obeyed by the proba- 
bilities (|55p and (|56p . since the scattering matrix is uni- 
tary; Eqs. ^ for the charge conservation are not obeyed. 
Following Refs.lolandllll. current conservation is now im- 
posed on Eqs. (O, leading to the determination of the 
chemical potential on the superconductor. This leads to 
a linear relation between II and the chemical potential 
difference fJ-L — fJ-R.i which is identified as the conductance. 
Denoting the latter by G^^, one has 



where^ 



(jr 



9i 



9ll9rr 9lr9rl 
9ll + 9rr + 9lr + 9rl 



2e 



x(*. 



sm' + \s., 



he 1 2 



(57) 



(58) 



Here, i and j refer to the two sides of the junction, say 
left and right, and the superscripts ee or he refer to the 
particular process. Thus for example, the 11 element of 
the matrix in Eq. ([52|) is 5£'^, while the 41 element is 

che 
'^RL- 

We compare the outcome of Eq. ([57)1 with the conduc- 
tance found from the Kubo formula in Fig. |4l There, the 
conductances are plotted for four values of the interface 
transmission, an almost perfect one, T = 0.96, (the up- 
most pair of curves), T = 0.8, (the second pair of curves 
from above), T = 0.5 and T = 0.31, (the low-lying two 
pairs of curves). In each case, the result of Eq. ([57|) 
is the dashed line. There are three interesting features 
of this comparison. Firstly, the conductance found from 
Eq. (|57[) is always smaller than the one found from the 
Kubo formula, Eq. (fT7)) . The difference between the two 
results decreases with increasing barrier (decreasing T) 



and seems to vanish in the limit T ^ 0. Thus, while for a 
normal conductor the ring geometry and the simple two- 
terminal configuration produce identical results for the 
conductance, this is unfortunately no longer the case for 
the NSN junction (NS for the ring geometry). The sec- 
ond interesting feature concerns the slopes of the curves 
in Fig. m when the barrier transmission, T is varied. 
While the conductance computed from the Kubo formula 
shows a crossover of the slope, from being positive at high 
values ofT to being negative at lower values, the slope of 
the conductance found from Eq. (j57p seems to be always 
negative. A third important difference between the two 
approaches is that for large ds/£, and not-too-small T, 
the Kubo result becomes larger than 1e} jh (tending to 
4e2//i in the limit 4/^ -> oo andT -> 1), while Eq. ((57)) 
actually tends to Gns = ^e^ jh and never yields the dou- 
bling of Gns due to the Andreev reflections. We blame 
these differences between the two approaches on the lack 
of conservation of charge in the BDG formulation. We 
believe that this deficiency is corrected by employing the 
Kubo formula for the ring geometry. 



G/(2e^/h) 
2.0 r 



0.5 



10 



20 



25 



30 



dM 



FIG. 4: Comparison between the conductance Eq. (I57p 
(dashed curves) and the conductance computed according to 
the Kubo formula, Eq. (|17|l (solid lines), as a function of 
ds/C for four values of the interface transmission. The latter 
conductance is always larger than the former. 

The difference between the two approaches becomes 
most marked in the limit of a nearly transparent barrier 
and a thick superconductor. This can be easily under- 
stood by noting that the addition of the two NS resis- 
tances is handled very differently by the two approaches. 
For the fully quantum case, adding two ideal conduc- 
tances (4e2/ft, for the NS case) gives just one ideal con- 
ductance. On the other, the scattering formalism, in the 
limit ds/C ~* o° and T — > 1 gives that both T (T') and 
7a (TJO vanish, and so does 7?. (7?.'), while 7?.a (T^a) 
tends to unity [see Eqs. (|55|l and ([56)) ]. As a result, 
9LL = 9RR = 4e2//i and gLR = 9r.l = 0, and Eq. ((57)) 
becomes exactly the classical addition of resistances, pro- 
ducing half the ideal quantum conductance of the pure 
NS junction. This is due to the fixing of the chemical 
potential on the S-section to conserve the current, as in 
the classical treatment. 
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In summary, we have shown that the presence of a su- 
perconducting segment in an otherwise normal system re- 
duces the overall conductance once the barriers between 
the superconducting and the normal parts become high 
enough. Thus, the superconducting segments may push 
the system towards the localizes insulating state. 

From the appearance of the plots presented in Figs. [H 
and [21 one may be tempted to say that the system expe- 
riences a metal-insulator quantum phase transition from 
a finite to a vanishing conductivity at large d^/Cj when 
T decreases (which can be inferred to as "disorder in- 
crease" ) . We refrain here from making such a statement 
and defer the discussion of such a quantum phase tran- 
sition in the thermodynamic limit for the composite NS 



system to future work. 
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